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ONE DIMENSIONAL DISSIPATIVE BOLTZMANN EQUATION: MEASURE 
SOLUTIONS, COOLING RATE AND SELF-SIMILAR PROFILE 


R. ALONSO, V. BAGLAND, Y. CHENG & B. LODS 


Abstract. This manuscript investigates the following aspects of the one dimensional dissipa¬ 
tive Boltzmann equation associated to variable hard-spheres kernel: (1) we show the optimal 
cooling rate of the model by a careful study of the system satisfied by the solution’s moments, (2) 
give existence and uniqueness of measure solutions, and (3) prove the existence of a non-trivial 
self-similar profile, i.e. homogeneous cooling state, after appropriate scaling of the equation. 
The latter issue is based on compactness tools in the set of Borel measures. More specifically, 
we apply a dynamical fixed point theorem on a suitable stable set, for the model dynamics, of 
Borel measures. 

Keywords. Boltzmann equation, self-similar solution, measure solutions, dynamical fixed 
point. 
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1. Introduction 

In this document we study a standard one-dimensional dissipative Boltzmann equation asso¬ 
ciated to "variable hard potentials" interaction kernels. The model is given by 

x) = Q(f, f)(t, x), (t, x) G [0, oo) X M, 

f(0,x) = fo(x), 

where the dissipative Boltzmann operator Q = is defined as 

Q{fJ){x)-= [ f{x-ay)f{x + by)\y\'^dy-f{x) I f {x + y) \y\'^dy . (1.2) 

M K. 

The parameters of the model satisfy 7 > 0, a G (0,1), 6 = 1 — a and will be fixed throughout 
the paper. Such a model can be seen as a generalisation of the one introduced by Ben-Na'im and 
Krapivsky [7] for 7 = 0 and happens to have many applications in physics, biology and economy, 
see for instance the process presented in [4, 7] with application to biology. 

The case 7 = 0- usually referred to as the Maxwellian interaction case - is by now well 
understood [17] and we will focus our efforts in extending several of the results known for that 
case (7 = 0) to the more general model (1.1). Let us recall that, generally speaking, the analysis 
of Boltzmann-like models with Maxwellian interaction essentially renders explicit formulas that 
allow for a very precise analysis [11, 13, 12, 29, 17, 7] because (1) moments solve closed ODEs 
and, (2) Fourier transform techniques are relatively simple to implement. In [7] a Brownian 
thermalization is added to the equation which permits a study of stationary solutions of (1.1) 
for 7 = 0. Here we are more interested in generalizing the works of [11, 13, 12, 29, 17] that 
deal with the so-called self-similar profile which in the particular case of Mawellian interactions 
is unique and explicit. Such a self-similar profile is the unique stationary solution of the self- 
similar equation associated to (1.1) and, by means of a suitable Fourier metric, it is possible 
to show (in the Maxwellian case) exponential convergence of the (time dependent) self-similar 
solution to this stationary profile, [13, 12, 17]. In particular, this means that solutions of (1.1) 
with 7 = 0 approach exponentially fast the "back rescaled" self-similar profile as t —)■ -|-oo. 

Self-similarity is a general feature of dissipative collision-like equations. Indeed, since the 
kinetic energy is continuously decreasing, solutions to (1.1) converge as t —)■ 00 towards a Dirac 
mass. As a consequence, one expects that a suitable time-velocity scale depending on the rate 
of dissipation of energy may render a better set up for the analysis. For this reason, we expect 
that several of the results occurring for Maxwellian interactions remain still valid for (1.1) with 
7 > 0. The organization of the document is as follows: we finish this introductory material with 
general set up of the problem, including notation, scaling, relevant comments and the statement 
of the main results. In Section 2 the Cauchy problem is studied. The framework will be the 
space of probability Borel measures. Such a framework is the natural one for equation (1.1) as 
it is for kinetic models in general. In one dimensional problems, however, we will discover that 
it is essential to work in this space in contrast to higher dimensional models, such as viscoelastic 
Boltzmann models in the plane or the space where one can avoid it and work in smaller spaces 
such as Lebesgue’s spaces [10, 32, 1, 27]. This last fact proves to be a major difficulty in the 
analysis of the model. The Cauchy problem is then based on a careful study of a priori estimates 
for the moments of solutions of equation (1.1) and standard fixed point theory. In Section 3 we 
find the optimal rate of dissipation (commonly referred as Haff’s law) which follows from a careful 
study of a lower bound for the moments using a technique introduced in [1, 2]. In Section 4 we 
prove the existence of a non-trivial self-similar profile which is based, again, on the theory of 
moments and the use of a novel dynamical fixed point result on a compact stable set of Borel 
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probability measures [5]. The key remaining argument is, then, to prove that the self-similar 
profile - which a priori is a measure - is actually an L^-function. Needless to say, such stable set 
is engineered out of the moment analysis of Sections 2 & 3. In Section 5, numerical simulations 
are presented that illustrate the previous quantitative study of equation (1.8) as well as some 
peculiar features of the self-similar profile G. The simulations are based upon a discontinuous 
Galerkin (DG) scheme. The paper ends with some perspectives and open problems related to 
(1.1), in particular, its link to a recent kinetic model for rods alignment [4, 8]. 


1.1. Self-similar equation and the long time asymptotic. The weak formulation of the 
collision operator Q reads 




j j f{x)f{y)\x-y\'^ (ij{ax + by) + 'ip{bx + ay)-'ip{x)-i;{y)^dxdy 


(1.3) 


for any suitable test function ip. In particular, plugging successively il>{x) = 1 and ip{x) = x 
into (1.3) shows that the above equation (1.1) conserves mass and momentum. Namely, for any 
reasonable solution f{t,x) to (1.1), one has 


f{t,x)dx 


fo{x)dx and 


/ xf{t,x)dx 

Jr 


xfo{x)dx Vt^O. (1.4) 


However, the second order moment is not conserved: indeed, plugging now ip{x) = |xp into (1.3) 
one sees that 

/ Q{fJ){x)\x\‘^dx = -ab[ /(x)/(2/) lx - yp+^dxdy 
Jr 4ir2 

since |ax -|- by\‘^ -I- \ay -|- 6xp — |xp — \y\‘^ = —2a b\x — yp for any (x, y) G M^. Therefore, for any 
nonnegative solution f{t,x) to (1.1), we get that the kinetic energy is non increasing 

^E{t) [ f{t,x)\x\‘^dx =-ab [ /(t, x)/(t, y) |x - yp+'^dxdy ^ 0 . (1.5) 

dr dt J]R2 

This is enough to prove that a non trivial stationary solution to problem (1.1) exists. Indeed, 
for any xq G M, the Dirac mass 5xq is a steady (measure) solution to (1.1). For this reason one 
expects the large-time behavior of the system to be described by self-similar solutions. In order 
to capture such a self-similar behavior, it is customary to introduce the rescaling 

y{t) 9{sit),0 = f{t:x), i = V{t)x (1.6) 

where V{t) and s{t) are strictly increasing functions of time satisfying H(0) = 1, s(0) = 0 and 
limi_^oo s{t) = oo. Under such a scaling, one computes the self-similar equation as 


x) = (V {t)g{s, 0 + ^ {t)s{t)dsgis, ^ + iV {t)d^g{s, 

(■5,6) + V{f)s{t)dsg{s 
while the interaction operator turns into 


s=s(t), ^=V(t)x 

s=s{t) , ^=V(t)x 


Q{fJ)(.t,x) = ^{t)Q{g,g){s{t),V{t)x). 


(1.7) 
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Consequently / = f{t,x) is a solution to (1.1) if and only if g = g{s,^) satisfies 
s{t)V'^{t)dsg{s,0 + = Q{g,g)is,(). 

Choosing 

1 1 

V{t) = (I + c^t)^ and s(t) = — log(l + C7t), c > 0 , 

C7 

it follows that g solves 

dsg{s,0 + cd^ {Cg{s,0) = Q{9,9 ){s,0 ■ ( 1 - 8 ) 

Thus, the argument of understanding the long time asymptotic of (1.1) is simple: if there exists 
a unique steady solution G to (1.8), then such a steady state G should attract any solution to 

( 1 . 8 ) and, back scaling to the original variables 

f{t, x) ~ C {t)G{V{t)x) as t —)■ oo . 

in some suitable topology. We give in this paper a first step towards a satisfactory answer to 
this problem; more specifically, we address here two main questions: 

Question 1. Determine the optimal convergence rate of solutions to (1.1) towards the Dirac 
mass centred in the center of mass xq := f^xfo{x)dx. The determination of this 
optimal convergence rate is achieved by identifying the optimal rate of convergence 
of the moments mk{t) of f{t,x) defined as 

RT'fc(i) := / \x— xq\^ f{t,x)dx k^O. 

Jr 

Question 2. Prove the existence of a “physical” steady solution G £ -^max (7 2 )(^) (^- 8 )) that 

is a function satisfying 

c^(^G(0) = Q(G,G)(e), eeK, (1.9) 

in a weak sense (where c > 0 is arbitrary and, for simplicity, can be chosen as c = 1 ). 
Note that equation (1.9) has at least two solutions for any 7 > 0, the trivial one 
and the Dirac measure at zero. None of them is a relevant steady solution since 
both have energy zero which is a feature not satisfied by the dynamical evolution 
of (1.8) (provided the initial measure is neither the trivial measure nor the Dirac 
measure). 

Similar questions have already been addressed for the 3-D Boltzmann equation for granular gases 
with different type of forcing terms [27, 23, 9]. For the inelastic Boltzmann equation in the 
answer to Question 1 is known as Haff’s law proven in [27, 1, 2] for the interesting case of 
hard-spheres interactions (essentially the case 7 = 1 ). The method we adopt here is inspired 
by the last two references since it appears to be the most natural to the equation. Concerning 
Question 2, the existence and uniqueness (the latter in a weak inelastic regime) of solutions to 

(1.9) has been established rigorously for hard-spheres interactions in [27, 28]. In the references 
[27, 23, 9], also [5, 20] in the context of coagulation problems, the strategy to prove the existence 
of solutions to problem similar to (1.9) is achieved through the careful study of the associated 
evolution equation (equation ( 1 . 8 ) in our context) and an application of the following dynamic 
version of Tykhonov fixed point theorem (see [5, Appendix A] for a proof): 
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Theorem 1.1 {Dynamic fixed point theorem). Let y be a loeally convex topological vector 
space and Z a nonempty convex and compact subset of y. If {J-t)t>o is a continuous semi-group 
on Z such that Z is invariant under the action of IFt (that is J-tZ G Z for any z G Z and t > 0), 
then, there exists Zq ^ Z which is stationary under the action of J-t (that is TtZo = Zq for any 
t > 0). 

In the aforementioned references, the natural approach consists in applying Theorem 1.1 to 
y = endowed with its weak topology and consider for the subset Z a convex set which includes 
an upper bound for some of the moments and some L^-norm, with p > 1, which yield the desired 
compactness in y. As a consequence, with such approach a crucial point in the analysis is to 
determine uniform L^-norm bounds for the self-similar evolution problem. This last particular 
issue, if true, appear to be quite difficult to prove in the model ( 1 . 8 ), mainly because the lack 
of angular averaging in the 1-D interaction operator Q as opposed to higher dimensional inter¬ 
action operators. This problem is reminiscent of related 1-D interaction operators associated to 
coagulation-fragmentation problems, for instance in Smoluchowski equation, for which propaga¬ 
tion of IP norms is hard to establish, see [24] for details. Having this in mind, it appears to us 
more natural to work with measure solutions and considering then T as a suitable space of real 
Borel measures endowed with the weak-* topology for which the compactness will be easier to 
establish. Of course, the main difficulty will then be to determine that the fixed point provided 
by Theorem 1.1 is not the Dirac measure at zero (the trivial solution is easily discarded by mass 
conservation). In fact, we will prove that this steady state is a function (see Theorem 1.5). 
Let us introduce some notations before entering in more details. 

1.2. Notations. Let us introduce the set A4s(M) as the Banach space of real Borel measures on 
M with finite total variation of order s endowed with the norm || • defined as 

||;u||s := / (x)® |;U |(dx) < oo , with (x) := (l-|-|xp) ^ Vx G M , 

Jr 

where the positive Borel measure \y\ is the total variation of p. We also set 

A4+(M) = {p e J^s(^); pP 0}, 

and denote by 'P(M) the set of probability measures over M. For any k y 0, define the set 

'Pfc(l^) = I M G V{R) ; / |x|^;u(dx) < oo 

I Jr 

For any p G 'Pfc(M) and any 0 ^ p ^ /c, let us introduce the p-moment Mp {p) := |x|^’p(dx). 

If p is absolutely continuous with respect to the Lebesgue measure with density /, i.e. p(dx) = 
f(x)dx, we simply denote Mp{f) = Mp{p) for any p P 0. We also define, for any k y 1, the set 

= \ h ^ ; [ xp{dx) = 0 

I Jr 

In the same way, we set Lj^(R) = L^(M) n A4 a:( 1^)) for any k y 0. Moreover, we introduce the 
set L^g {s y 0) of locally bounded Borel functions (p such that 

:= sup |(p(x) I (x)“® < oo . 
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For any p ^ 1 and £ 'Pp(M), we recall the definition of the Wasserstein distance of order p, 
Wp{p, u) between p and u by 

Wp{p,u)=( irif [ |x - y|^’ 7 r(dx,dy)^ , 

VTi-eni/i,!/) J^2 J 

where n(|U, u) denotes the set of all joint probability measnres vr on whose marginals are p 
and u. For the pecnliar case p = 1, we shall address the first order Wasserstein distanee Wi{p, v) 
as Kantorovieh-Rubinstein distanee, denoted c^kR) be- dKR{p,i^) = Wi{p,iy). We refer to [33, 
Section 7] and [34, Chapter 6 ] for more details on Wasserstein distances. The Kantorovich- 
Rnbinstein dnality asserts that 


dKR{p, = snp 

veLipiO 


<p{x){p - n){dx) 


where Lipj 


denotes the set of Lipschitz fnnctions p snch that 

\p{x) - <p{y)\ 


= sup 

x^y 


\x - y\ 


^ 1 


For a given T > 0 and a given /c ^ 0, we shall indicate as Cweak([0, T], 'Pfc(M)) the set of continuous 
mappings from [0,T] to 'Pfc(l^) where the latter is endowed with the weak-* topology. 


1.3. Collision operator and definition of measure solutions. We extend the definition 
(1.3) to nonnegative Borel measures; namely, given Pjf £ Af;Jr(M), let 

{Q{p,u)-, ip) :=]- [ \x-y[' l^ip{x,y)p{dx)v{dy) (1.10) 

4 Jr2 


for any test function tp £C (M) H (M) where 

A(p{x,y) :=ip{ax + by) + ip{bx + ay) - (p(x) -(p(y), (x,y) G . 

A natural definition of measure solutions to (1.1) is the following, (see [25]). 


Definition 1.2. Let 7 > 0, T = max( 7 , 2 ), po £ A4y(^) (^t )«70 C A4y(J^) given. We 

say that {pt)t^o is a measure weak solution to ( 1 . 1 ) associated to the initial datum po if it satisfies 

(1) supt^o \\Lt\\r ^ lllioIlT, 



= / Ro{dx) 
Jr 


and 


/ xpt{dx) = / xpo{dx) Vt > 0 . 

M. M 


( 2 ) for any test-funetion (p £ C;,(M) := C(M) H L°°(M), the following hold 

i) the mapping 1 1 —)■ {Q{pt, Pt) ] p) belongs to C([ 0 ,oo)) C ([ 0 ,oo)), 
a) for any t 0 it holds 


(p{x)pt{dx)= / (p{x)po{dx)+ / {Q{pr,PT) ; p)'iT. 
: Jr Jo 


Notice that if pt £ ^4^(1^) for ^-riy t 0, then 


( 1 . 11 ) 


( 1 . 12 ) 



^ |A(p(x,p)| pt{dx)ptidy) ^ 4||(p||oo||pi||^ 


< 00 


for any p £ Cb(K) and t 0. This shows that {Q{pt, pt)fi) is well-defined for any tfid and any 
ip £ Cb(M). Similarly, the notion of measure solution to (1.9) is given in the following statement. 
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Definition 1.3. A measure G ^max (7 2 )(^) ^ solution to (1.9) if 

- [ = A / \C-Vr{<P{af + b'n) + (f{ar] + b^)-(j){f)-(j){rj))n{df)fi{dr]) (1.13) 

Jr ^ Jr2 

for any (j) G Cb(K) where (f' stands for the derivative of (j). 

Notice that, by assuming \\n\\o = 1, we naturally discard the trivial solution G = 0 to (1.9). 


1.4. Strategy and main results. Fix 7 > 0, T = max( 7 ,2). Thanks to the conservative 
properties ( 1 . 11 ), we shall assume in the sequel, and without any loss of generality, that the 
initial datum ^0 £ A4y(M) is such that 


//o(da;) = 1 and / xfio{dx) = 0 , 

; Jr 


i.e. /io £ 'Py(^)- This implies that any weak measure solution to (1.1) associated to ^uq 

is such that 


^jGiP^(M) Vt^O. 


Theorem 1.4. Let fj,Q G 'Py(^) be given initial datum. Then, there exists a measure weak 
solution {pt)t^o to (1.1) assoeiated to fiQ in the sense of Definition 1.2. Moreover, 



^o(dx) = 1, 


/ x^i(dx) = / xfio{dx) = 0, 
«/ M M. 


and 



^htidx) 




x|^/ro(dx) 


Vt ^ 0. 


Moreover, if fio £ nA: 7 oT’fc(^) then C nfc 7 oT’fc(I^)- U additionally there exists e > 0 sueh 
that 

/ exp(e|xP)/io(dx) < 00 , (1-14) 

Jr 

then, sueh a measure weak solution is unique. Furthermore, if no is absolutely eontinuous with 
respeet to the Lebesgue measure, i.e. /Uo(dx) = /o(x)dx with /o G Ly(M) then pLt is absolutely 
eontinuous with respeet to the Lebesgue measure for any t ^ 0. That is, there exists {ft)t^o C 
LyfR) sueh that p.t{dx) = ft{x)dx for any t ^ 0. 


We prove Theorem 1.4 following a strategy introduced in [25]-[22] for the case of Boltzmann 
equation with hard potentials (with or without cut-off). The program consists essentially in the 
following steps. (1) Establish a priori estimates for measure weak solutions to (1.1) concern¬ 
ing the creation and propagation of algebraic moments, ( 2 ) construct measure weak solutions 
to (1.1) by approximation of L^-solutions. Step (1) helps proving that such approximating se¬ 
quence converges in the weak-* topology. Finally, (3) for the uniqueness of measure weak solution 
we adopt a strategy developed in [22] and based upon suitable Log-Lipschitz estimates for the 
Kantorovich-Rubinstein distance between two solutions of (1.1). 


As far as Question 1 is concerned, we establish the optimal decay of the moments of the 
solutions to (1.1) by a suitable comparison of ODEs. Such techniques are natural and have been 
applied to the study of Haff’s law for 3-D granular Boltzmann equation in [1]. The main difficulty 
is to provide optimal lower bounds for the moments, see Propositions 3.3 & 3.4. Essentially, we 
obtain that 
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see Theorem 3.5 for a more precise statement. Such a decay immediately translates into conver¬ 
gence of towards ho in the Wasserstein topology. 

Regarding Question 2, once the hxed point Theorem 1.1 is at hand, the key step is to engineer 
a suitable stable compact set Z. Compactness is easily achieved in T = ^max( 7 , 2 )(^) endowed 
with weak-A topolgy; only uniformly boundedness of some moments suffices. However, y must 
overrule the possibility that the hxed point will be a plain Dirac mass located at zero. This 
is closely related to the sharp lower bound found for the moments in Question 1. In such a 
situation, a serie of simple observations on the regularity of the solution to (1.13) proves that 
such a steady state is actually a L^-function. Namely, one of the most important step in our 
strategy is the following observation: 

Theorem 1.5. Any steady measure solution fi £ ^max (7 2 )(^) (1-9) sueh that 

1^7 ;= / |^r/^(dO > 0 (1.15) 

Jr 

is absolutely continuous with respect to the Lebesgue measure over M, i.e. there exists some 
nonnegative G £ .^"max (7 2 ) (^) that 

tx{do = G{m. 

In other words, any solution to (1.9) lying in ^max (7 2 )(^) different from a Dirac mass must 
be a regular measure. This leads to our main result: 

Theorem 1.6. For any 7 > 0, there exists G £ .^"max (7 2 ) which is a steady solution to (1.9) 
in the weak sense. 


2. Cauchy theory 

We are hrst concerned with the Cauchy theory for problem (1.1) and we begin with studying 
a priori estimates for weak measure solutions to (1.1). Let us hx 7 > 0 and set T = max( 7 , 2). 


2.1. A priori estimates on moments. We hrst state the following general properties of weak 
measure solutions to ( 1 . 1 ) which have sufficiently enough bounded moments. 

Proposition 2.1. Let po £ AIy(I^) let be any weak measure solution to (1.1) 

associated to //q. Given k y 0, assume there exists p > k -£ 'j such that 


sup ||//i||p<oo, y T > 6 > 0. (2.1) 

S<t<T 

Then, the following hold: 

(1) For any ip £ L“^(M) nC(M), the mapping t ^ 0 1 —)■ {Q{pt, Pt) ; t) continuous in (0, 00 ). 

(2) For any p £ L5)().(M) C C(M) it holds that 

d f 

— p{x)pt{dx) = {Q{pt, Ft)-, t) ytyO. ( 2 . 2 ) 

at Jr 

The proof of Proposition 2.1 will need the following preliminary lemma (see [25, Proposition 
2 . 2 ] for a complete proof). 
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Lemma 2.2. Let he a sequenee from eonverges weakly-k to some ^ G 

We assume that for some p > 0 it holds 


'ijj(x y\ 

Then, for any f; G C(M^) satisfying lim — - ’ = 0, one has 

IxI+Ij/K+oo {x)p + {y)P 

lim / 'ip{x,y) p^-idx) p'^idy) = i>{x,y) p{dx) p{dy). 

n^+oo J^2 J^2 

Proof of Proposition 2.1. Let /c ^ 0 and p G n C(]R) be given. Choose x £ such 

that llxiloo ^ 1 with x{x) = 0 if |x| ^ 2 and xi^) = 1 for 1^1 ^ 1 ^-rid set for any n G N*, 
Tn{x) = p{x)x (^)- It follows that pn £ Cc(M) C Cf,(M) for any n with Pn{x) —)• p{x) for any 
X G M as n —)■ oo. Consequently, Apn{x,y) —)■ Ap{x,y) for any {x,y) G as n —)■ oo. Now, 
since pn G C;,(M) for any n G N, one deduces from (1.12) that 

f Pn{x)pt2{dx) = f Pn{x)pti{dx) + f {Q{pr, Pt)T n) dx \/ t2 > h > 0 , Vn^l. 

Jm. Jm. Jti 

Notice that ^ < oo. Thus, there is C > 0 such that for any n G N*, any 

{x,y) G M^, 

|Av9(x, 2/)| ^ C (^(x)^ + and \Apn{x,y)\ ^ C (^{x)’^ + {y)'^'^ , (2.3) 

where the constant C depends only on p (and a). Using the dominated convergence theorem 
together with (2.1), one deduces that 

rL P 2 r 

/ {QiprPr); Tn)dT = - dx \x - yl"^Apn{x, y)Pr{dx)Pr{dy) 

Jt\ ^ Jti J 

rt 2 r rt 2 

—> - / dx \x - ypAp{x,y)pr{dx)pridy) = / {Q{pr, Pr); p) dx , 
n^cx^ 2 J^2 Jt 


SO, the identity 


p{x)pt2{dx)= f p{x)pt^{dx)+ f {Q{pr,PT); T)dx V t2 > U > 0 , ( 2 . 4 ) 

Jm. Jti 


holds true. It follows from (2.3) that 

{Q{Pt, Pr) ;t) J^Jx- yP (^(x)^ + {y)^'j pr{dx)pr{dy) ^ 2C\\pr\\l+^ Vr > 0 . 
Combining this with (2.4), one sees that, for any T > <5 > 0, 

/ p{x)pt 2 {dx)- / p{x)pti{dx) ^2C'|f2-U| sup \\pT\\l+.y y h, t 2 e [5,T] . 

Jr Jr s^t^t 

In particular, under assumption (2.1), the mapping 1 1 —)■ / p{x)pt{dx) is continuous over (0, co). 

Jr 

This shows that, for any t > 0 and sequence (fn)n C [f/2,3t/2] with lim„t„ = t, the sequence 
{pt„)n converges weakly-* towards pt. In addition, since p G and p > k + x one concludes 

that 

\x - y\^Ap{x,y) 

hm ----- = 0 . 

|a;| + |j/Koo {x)P -h {y)P 
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Thus, it readily follows from Lemma 2.2 that 


lim 

n 



X - 2/PAv9(x,y)/rt„(dx)//t„(dy) 


X - A(^(x,y)//t(dx)^t(dy). (2.5) 


Henceforth, the mapping t i—)■ is continuous over (0,oo) proving point (1). Point 

(2) follows then directly from (2.4). □ 

Moments estimates of the collision operator are given by the following: 

Proposition 2.3. Let /i G 'Pfc+..y(M) with k ^ 2. Then 

(2.6a) 

In partieular, 

{Qip, /i); I • I") ^ (l - a" - 6") Mfc (/i)i+i , (2.6b) 

and, if k > 2 

{Q{t,p) ; I • I") ^ - 2 (l - a" - b^) M 2 (m)“* Mfc(^)^+^. (2.6c) 

Proof. We apply the weak form (1.10) to ip{x) = \x\^. Using the elementary inequality 

|x|^ + \y\^ - \ax + by\^ - \ay + 6x|^ ^ (^1 - - b^^ |x - y\^ (2.7) 

valid for any (x,y) G and k ^ 2 (with equality sign whenever k = 2), and noticing that 
1 — — 6^ is nonnegative for any k ^ 2 and any a G (0,1) we have 

{Qih, /^); I • I") ^ (l - a" - b^) \x - yr+V(dx)^(d 2 /). (2.8) 


Since /a G P^_|_.^(M) and the mapping [0, oo) 9 r i—)■ is convex, one deduces from Jensen’s 

inequality that 


/ k - y\'^^^Tidx)fi{dy) 
Jr2 




x|^^'^;u(dx) , 


from which inequality (2.8) yields (2.6a). Using again Jensen’s inequality we get also that 
Mfc_|_.y(^) ^ which proves (2.6b). Finally, according to Holder inequality 


Mk+^{p.) ^ Mfc(/i)^+'=-2 M2iy) '“-2 yk>2, 


and (2.6c) is deduced from (2.6a). 


□ 


One deduces from the above estimate suitable upper bounds for the moments of solution to 

( 1 . 1 ): 

Proposition 2.4. Let fj,Q G be a given initial datum and {iat)t^o be a measure weak 

solution to (1.1) assoeiated to fiQ. Then, the following holds 

(1) If iPt) t^o C Pp(K) with p > k + ^ for some k ^ 2, then, t ^ 0 i—)■ Mk{y.t) is deereasing 
with 

_ k 

Mk{^it)^Mk{^io)(l + ^{l-a^-b>^)Mk{^ioy^t) " Vt^O. (2.9) 

In partieular, sup^^Q < oo. 
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(2) If for some k > 2 and p > k + ^, \\l^t\\p < oo for any (5 > 0, then 

Mkipt) ^ Ck{'y)\\lJ,o\\T m.m{t Vt>0, (2.10) 

where Ck{'y) > 0 is a positive constant depending only on k > 2, ^ and a. 

Proof. For the proof of (2.9), apply (2.2) with p{x) = \x\^ and note that using (2.6b) 

^Mfc(/it) ^ -1 (l - - 6 "=) Mk (^i)^+fc Vt ^ 0 . 

This directly implies the point (1) of Proposition 2.4. To prove estimate (2.10) for short time 
observe that applying ( 2 . 2 ) to (p{x) = \x\^ and using ( 2 . 6 c), one gets 

^Mfc(^i) ^ -1 (^1 - - 6 ^) M 2 (//t)-s^ Mfc(Mi)i+E^ V t > 0 , 

and, since by definition of measure weak solution it holds M 2 {pt) ^ IImoIIt for any t ^ 0 , one 
deduces that 

jMkipt) llliollx^ Mfc(/it)'+^ V t > 0 . 

This inequality leads to estimate (2.10) with constant Ckij) = ( 2 {k- 2 ) (l ~ j ■ The 

long time decay follows applying (2.9) for t ^ 6. □ 

Remark 2.5. It is not difficult to prove that the conclusions of the previous two propositions 
hold for general measure po G A4y(IP) with 

llirollo = £» / 0 and / x^o(da;) = 0 . 

Jr 

In such a case, the constant Ck{'j) depends also continuously on g. 


Introduce the class V, 


exp, 7 


of probability measures with exponential tails of order 7 , 


Ve 


exp, 7 v^^i = G 'P(M) ; Be > 0 such that J exp(e|xP);u(dx) < 00 j . (2-11) 

We have the following: 

Theorem 2.6. Let po G 'Py(^) initial datum and {pt)t bo o measure weak solution to 

(1.1) associated to pQ. If po G 'Pexp, 7 (l^) ond if {pt)t C P|^^q'P)?(M), then, there exists a > 0 and 
C > 0 such that 


sup / exp(a|xP)/it(dx) ^ C. 
*70 Jr 


( 2 . 12 ) 


Proof. Since po G Vexp,'y(^), there exists a > 0 and Cq > 0 such that 

/ exp(Q;|xP)/io(dx) ^ Cq . 

Jr 


Let us denote by po the integer such that 'ypofT and 7 (po — 1) < T. Thus, for 0 ^ p ^ pQ — 1 
and t f 0, 

M^p{pt) f llliillx ^ lll^ollx- 
For p ^ Po and t f 0, one deduces from (2.6a) that 

M^piyPt) f M^pi^po'). 
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Thus, for any f ^ 0 and any n > po 

n p Po -1 p n p 

^ II/xoIIt ^ ||//o||Texp(a) tCq. 

p=0 p=0 p=po 

Letting n —)• +oo we get (2.12). □ 

2.2. Cauchy theory. The main ingredients of the proof of Theorem 1.4 are Propositions 2.7 
and 2.8. The details of the proof can be found in Appendix A. Namely, by studying hrst the 
Cauchy problem for initial data and then introducing suitable approximation we can construct 
weak measure solutions to (1.1) leading the following existence result: 

Proposition 2.7. For any po £ /^o 7^ there exists a measure weak solution C 

'Py(^) associated to po and such that 

sup ll^tllx ^ IImoIIt and sup ||;Ut < oo Vto>0, s>T. 

Moreover, if po £ {pt)t;»o C 

To achieve the proof of Theorem 1.4, it remains only to prove the uniqueness of the solution. 
It seems difficult here to adapt the strategy of [25] for which the conservation of energy played a 
crucial role. For this reason, we rather follow the approach of [22] which requires the exponential 
tail estimate (1.14). The main step towards uniqueness is the following Log-Lipschitz estimate 
for the Kantorovich-Rubinstein distance. 


Proposition 2.8. Let po and uq be two probability measures in 'P 2 (^) satisfying (1.14) i.e., 
there exists e > 0 such that 


/ exp (e |xP) ^o(dx) + 

Jm 


exp (e |xP) i^o(da:) < oo . 


Let then {pt)t^o and be two measure weak solutions to (1.1) associated respectively to the 

initial data po and vq. There exists > 0 such that, for any T > 0, 

dKRiTt,t^t) ^ dKR{po,i^o) (2.13) 

+ KsCT{e)J (iKR(/^s,i^s)(l + I logdKR(lis,z^s)|)ds VtG[0,r], 

where C'r(e) = sup / exp (e |xp) {pt + vt) (dx) < oo. 
te[o,r] Jr 


Proof of Theorem 1.4- Given Proposition 2.7, to prove Theorem 1.4 it suffices to show the unique¬ 
ness of measure weak solutions to (1.1). Let po be a probability measure in satisfying 

(1.14) and let {pt)t^o and t>e two measure weak solutions to (1.1) associated to pQ. From 

Proposition 2.8, given T > 0 there exists a hnite positive constant Kt such that 

rt 

is Vf G [0, T] , 


dKR{ht,vt) ^ Kt / $(dKR(/ns, J^s))ds 
Jo 


with $(r) = r (1 -|- I logr|) for any r > 0. Since $ satishes the so-called Osgood condition 

dr 

/ ^ = oo, (2.14) 

Jo ^{r) 

then, a nonlinear version of Gronwall Lemma (see for instance [6, Lemma 3.4, p. 125]) asserts 
that dKR{pt,i^t) = 0 for any t £ [0,r]. Since T > 0 is arbitrary, this proves the uniqueness. □ 
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The existence and uniqueness of weak measure solution to (1.1) allows to define a semiflow 
(^t) on 'Pexp, 7 (lK), (recall (2.11)). Namely, Theorem 1.4 together with Theorem 2.6 assert 
that, for any fiQ G 'Pexp, 7 (l^) H 'P^(M), there exists a unique weak measure solution (^t)t^o to 
(1.1) with fj,t G 'Pexp, 7 (l^) for any t ^ 0 and we shall denote 

l^t-=St{no) Vt^O. 

Then, the semiflow St is a well-defined nonlinear mapping from 'Pexp, 7 (I^) H 'P^(M) into itself. 
Moreover, by definition of weak solution, the mapping 1 1 —5t(^o) belongs to Cweak([0, oo), ^ 2 ( 1 ^)). 
One has the following weak continuity result for the semiflow. 

Proposition 2.9. The semiflow {St)t^o is weakly continuous on 'Pi(M) in the following sense. 
Let {ian)n G 'Py(^) ® sequence such that there exists e > 0 satisfying 

sup / exp(e|xP)/rn(dx) < 00 . (2.15) 

nSN JR 

If {nn)n converges to some fi G 'Py(^) t/ie weak-k topology, then for any t ^ 0, 

StilXn) —^ <Stih) t/ie weak-k topology as n —)• 00 . 


Proof. The proof is based upon the stability result established in Proposition 2.8. Namely, 
because {p,n)n converges in the weak-* topology to ^ G 'Py(^)) deduces from (2.15) that 

/ exp(e|xP)^(dx) < 00 , 

Jr 

which means that all p,n and fx share the same exponential tail estimate with some common 
e > 0. Then, for any T > 0 one deduces from (2.13) that 

^^KR(‘5t(//n) )‘5t(^)) ^ dKR(/^n i h) 

+ K{£)CT{e)[ ^{dKRiSsifXn), ■Ss{fi))ds Vt G [0,T], n G N, 

Jo 

for some universal positive constant K(e) and 

Crie) = sup sup / exp(e|xP) {St{pLn) + Stin)) (dx) < 00 
neNte[o,r] Jr 

according to Theorem 2.6. Here above, $(r) = r (1-|- |logr|) satisfies the Osgood condition 
(2.14), thus, using again a nonlinear version of Gronwall Lemma [6, Lemma 3.4] we deduce from 
this estimate that 


{Pn , h) ) - ^(^KR {•StiiXn) , ‘5t(^))) ^ iL(e) Ct{s) t Vt G [0,r], n G N 


where 


T'(x) = 


dr 

$(r) 


Vx ^ 0. 


Taking now the limit n —)■ 00 , since c^kr metrizes the weak-* topology of 'Pi(M) it follows that 
lim„^oo dKR{p.n, h) — 0- Furthermore, recalling that ^'(O) = 00 one concludes that 


lim ^(dKRiStiUn), Stih))) = 00 . 

n^oo ^ ^ ^ 

That is, lim„^.oo dKR(‘5t(^n), •St{P‘)) = 0 which proves the result. 


□ 
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3. Optimal decay of the moments 


We now prove that the upper bounds obtained for the moments of solutions to (1.1) in Propo¬ 
sition 2.4 are actually optimal. 

3.1. Lower bounds for moments. We begin with the case 7 G (0,1]: 

Proposition 3.1. Fix 7 G (0,1] and let fiQ G initial datum and be a measure 

weak solution to (1.1) assoeiated to Then, there exists JC^y > 0 depending only on a and 7 
sueh that 

JOw) > , , Vi>0. (3,1) 


Thus, 


1 -I- IC,yMy{HQ)t 
M2{fit) > ( 1 + IC^ M,y{p,o)t 


Vt ^ 0 . 


(3.2) 


Proof. For 7 G (0,1] the following elementary inequality holds 

B^{x, y) = (^ \ax + \ay + bx\'^ - \x\'^ - \y\'^'^ |x - ^ -C^\x\'^ |yP Vx, y G M, (3.3) 

for some positive constant C.y > t) explicit depending only on a and 7 . Using (2.2) and the weak 
form (1.10) with (p{x) = \x\^ G L^ 2 i^) nC(M), we get that 

d 1 /* C f 

—= - / B.y{x,y)nt{dx)ixt{dy) ^\x\'^ \y\^yt{dx)y.t{dy) 
dt 2 J^2 2 J^2 

from which it follows that 

VOO. 

This inequality, after integration, yields (3.1) with Inequality (3.2) follows from the 

fact that according to Holder inequality 

^ M2{ixt)2 WixtWQ^ , 

while ll^illo = ll^ollo = 1 for any t ^ 0. □ 

The remaining case 7 > 1 is more involved. In this case, inequality (3.3) no longer holds and 
we need the following result. 

Lemma 3.2. Fix 7 > 1. For any p > 1 and y G Fp_|_,.^(M) it holds that 

{Q{t, m) ; I • D ^ - (1 - l3p{a))Mp+^{n)+ np{fx), (3.4) 

where f32p{a) = max{a^“^, and 

[^] 

^p(f) ~ Ppiof) ^ ^ ( A; ) fc(F) + fo^fc(M) fo^p-fc+7(/^)^ • 

k=l 


Moreover, for k G (0,1] 

-{Qip , f) ; I • 1^) ^ C^,k M^{y) Mkipt). 

The eonstant C^^k depends only on 7 and k. 


(3.5) 
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Proof. Let us begin with inequality (3.4). We first notice the following elementary inequality, 
valid for any p > 0, 

\ax + by\P + \ay + bx\P ^ f3p{a) + y^)^ ^ /?p(a) (|x| + \y\f V(x, y) G , (3.6) 

where k i—)• /3k{a) is decreasing with / 32 (a) = 1 and lim^^oo /3fc(a) = 0 (recall that max{a, 6 } < 1). 
We use then the following useful result given in [14, Lemma 2] for estimation of the binomial for 
fractional powers. For any p ^ 1, li kp denotes the integer part of the following inequality 
holds for any u,v G M+, 


{u + vY-uP-yP [u^yP-^ + yP-^y^ , 


(3.7) 


k=l 


where the binomial coefficients are defined as 


/ p \ r p{p-i)-{p-k+i) ^ ^ ^ 

y 1 for/c = 0 . 

Therefore, for any p > 1, one gets 

|ax + hy\P + |ay + bx\P - \x\p - \y\P ^ -(l - /3p(a)) {\x\p + |y|^) 


+ft(<.)^(y)(ix|'=i!,r‘+iir'=i9i'=). 

k=l 

Therefore, using (1.10) we get 

(Q(/^, m); I • D ^f (kr + l2/r)k-y|V(dx)p(dy) 

4R2 

ft 

+ \^p{a) ^ ^ - yr + \x\P~'^ |y|^)p(dx)p(dy). 

Furthermore, since 7 > 1 using Jensen’s inequality and the fact that p G 'P®(M), one obtains 

(kr + lyD \x - y\^p{dx)p{dy) ^ 2 Mp+^{p). 

Consequently, 

{QiP, /^); I - D ^ -( 1 -^p(a))Mp+^(/r) 

^P ft 

+ I /3p(a) E(l)£jx-yrHxl^ |yr^ + \x\P-^ \y\Y p{dx) p{dy). 

For the remainder term, one uses the inequality \x — y|'>' ^ 2'>'“^(|x|'^ + ly]"'') to get that 

\x - yr {\x\^\y\^~^ + \x\^~’^ \y\Yp{dx)p{dy) 

^ 2 ^(Mfc+^(p) Mp_k{p) + Mfc(y) Mp_k+Yp)) > 

which proves (3.4). The proof of (3.5) relies on inequality (3.3). Indeed, for any p and any 
k G (0,1) one has 

{QiP ’ ’ I ■ 1^) = ^ ( I®* + “ |yl^) I* “ 2/r /a(dx)y(dy). 
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Now, by virtue of (3.3) with k G (0,1) instead of 7 

\ax + by\'" + \ay + - |x|^ - \y\^'^\x - y\^ ^ -Ck\x\^ \y\^ \/x,y gR, 

for some positive constant > 0 explicit and depending only on a and k. Therefore, 

-{Qiy , ; I • I") ^ ^ / \x\^ |y|" k - yr^fiidx^dy) . 

Using that \x — y\'^~^ ^ max{l, 2 '!'“^“^} {\x\'^~^ + \y\'^~’^) we get 

-{Qifi , ;u); I • I") ^ ^ max{l, 2^-^-^} |x|" \y\^ (ix^-^ + {y^^) ^(dx)^(dy) 

= max{l, 2 '^“*'“^} CkM^{n)Mk{n) 

which gives the proof with = max{l, 2 '*'“^“^} Ck- □ 

As a consequence, the following proposition. 

Proposition 3.3. For any 7 > 1 and s G (0,1] it follows that 

Ms+^iHt) ^ ^ ^ 0, (3.8) 

with As any constant such that 

/I ^ //A Ms+'riho) \ /Q 

As ^ max \ C'y^sy i_i _2 i ’ (^■^) 

[ Ms^^yo)) 

and C^^s is a constant depending only on 7 and s. 

Proof. Define X{t) ;= Ms+'^{p,t)—AMs{ytY~^^ where the constant A will be conveniently chosen 
later on. Since 7 + s > 1 and s G (0,1] we can use Lemma 3.2 to conclude that 

^X(t) = ^^Ms+,{yt) - ^ (1 + ^) Ms{ytfA^Ms{^it) 

^ — (1 — I3s+'y{a)) Ms+2'iiht) + 'F's+'yipt) + X ■ (3.10) 

Let us observe that for such a choice of s and 7 , one has < s +7 and for 1 ^ A: ^ , 

a simple use of Holder and Young’s inequalities leads to, for any e > 0, 

fc+T' s+7' —fe s+27 

Mk+yint) Ms+y-kiht) ^ Ms+2^{yt)‘^ Ms+^{yt) ^ eMs+2'y{Ft) + Af^+^(/rt)^^ 

with Kg^'y^e > 0 a constant depending only on s, 7 and e. Similar interpolation holds for the 
terms of the form Mk{yt) Ms+ 2 'y-kiFt) appearing in TZs+'yiht)- As a consequence, for any e > 0 , 
we have 

■S+27 

^ ^-^5,7-^s+27(/^0 H“ -^ 5 , 7 -^s, 7 ,e-^s+ 7 (Mi) ? (3.11) 

where > 0 is a constant depending only on s and 7 . Choose e = 6 ( 7 ) > 0 such that 

2 Xg^^ e ^ 1 /? 5 -|-^(tt) 

we obtain from estimates (3.10) and (3.11) that there are constants Ky = Ki{s,’j) and K 2 = 
7^2 (- 5 , 7 ) such that 

—X[t) ^ —— (1 — /3s_|_^(a)) Mg^2'y{Ft) 

+ AKiMs{yt)^'y Mg+^{iat) 'i + K 2 Mg+^{iat), 


(3.12) 
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where we also used the interpolation 
As a final step, notice that 

Msint) ^ ^"+ 27 (^ 4 ) • (3-13) 

Therefore, including (3.13) in (3.12) one finally concludes that 

—X(t) ^ —— (1 — /3s+-y(a)) Ms+'^{nt) 

1 I 7 I s 1 s s+27 

+ -< -< + K 2 ■ (3.14) 

Now, choosing A > 0 such that X(0) < 0, if there exists to > 0 for which X{to) = 0, then 
estimate (3.14) implies that 

^A(to) ^ ■ (3.15) 

Then, choosing A = A( 7 ,s) sufficiently large such that the term in parenthesis in (3.15) is 
negative we conclude that X'{to) ^ 0. This shows that, for such a choice of A, X{t) ^0 for any 
f ^ 0. □ 


Using Proposition 3.3 the desired lower bound is obtained. 


Proposition 3.4. For any 7 > 1 and s G (0,1] one has 

m.(h) >-^ , 

{1 + C„yl. 7 Af.(,«)!()’ 

where depends only on 7 and s and Ag is given by (3.9). Moreover, it holds 

Mkint) ^ ^ V/c > 1, s G (0, 1] and t ^ 0 

(1 + C^,,A, "2M,(/io)?t)- 

Proof. Using Lemma 3.2 and inequality (3.8), since s G (0,1) one gets 


V /c > 1, s G (0,1] and t ^ Q. 


-—Msint) = I • r) ^ Msifit) M.y{iJ,t) 

^ C^,,M,(/it)^+7 Mg+.yiiat)^-^ ^ C.y,sAl~^ M,(^t)^+? , 
where, for the second estimate, we used the inequality 

Integration of this differential inequality leads to 

Msilit) > --T . 

(1 + C^,,A, - ^ Msiiao)^ t)- 
The second part of the result follows from 


Msipt) ^ MUpt) Moiptf—^ = MM Pit), 


(3.16) 


(3.17) 


valid for any k > 1. 
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We can now completely characterize the decay of any moments of the weak measnre solntion 
associated to fiQ. Namely, 

Theorem 3.5. Let k ^ 2 and hq G be given. Denote by o. measure weak solution 

to (1.1) assoeiated to the initial datum Then, 

(1) When 7 G (0,1], there exists some universal constant /C.y > 0 (not depending on hq) such 
that for any p ^ 7 

£ 

- Vt^O, (3.18a) 

^1 + /C 7 Afly(/i 0 ) t^ 

and for any p G [2,k] 

Mp {pit) ^ -(Mm) ---- V t ^ 0 . (3.18b) 

{l + ^{l-aP-hP)Mp{piQ)'i ty 

In particular, if piQ G Vexp,'r{^) then (3.18b) holds true for any p ^ 2. 

(2) When 7 > 1, for any p > I, 

-- jr^Mpifit) Vt^O, (3.19a) 

(i+7C'^,iA Mi(^o)^t)" 

and for any p G [2,k] 

Mp {lit) ^ -MM---- V t ^ 0 , (3.19b) 

{l + q-^{l-aP-hP)Mp{iio)'v ty 

where C.y^i and Ai are defined in Proposition 3.4- 

Proof. The npper bonnds in (3.18b) and (3.19b) have been established in Proposition 2.4. For 
the lower bonnd, whenever 7 G (0,1), one simply nses Jensen’s ineqnality to get 

Mp{iit) ^ M^{lit)^ \ftft0, p^7, 

and, then conclnde thanks to eqnation (3.1). For 7 > 1, the lower bonnd is jnst (3.17) with 
s = 1. □ 

Corollary 3.6. Fix 7 > 0 and let hq G 'P®(M) n'Pexp,7(i^) bs an initial datum. Then, for any 
p ^ 1, the unique measure weak solution {lit)t^o is converging as t —)■ 00 towards Sq in the weak-k 
topology o/'Pp(M) with the explicit rate 

Wp{iit,do)ix{l + Ct)-- as t —)• 00 , 
for some positive constant C depending on fiQ. 

Proof. The resnlt is a direct conseqnence of Theorem 3.5 since 

Wp{iit,5o)P= [ \x\Piit{dx) = Mp{iit) t^O, 

Jr 

and Wp metrizes 'Pp(M) (see [34, Theorem 6.9]). □ 
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3.2. Consequences on the rescaled problem. Given /Uq £ 'Px(^) satisfying (1.14) and let 
denotes the unique weak measure solution to (1.1) associated to /tq- Recall from Section 
1 the definitions of the rescaling functions 


V{t) ■.= {1 +and 


s(t) := — log(l+ C 7 t), (c> 0 ). 

C7 


For simplicity, in all the sequel, we shall assume c = 1. The inverse mappings are defined as 

exp (7 s) — 1 


V(s) = exp(s) t{s) = 


7 


i.e. 

s{t) = s t{s) = t and V{t{s))=V{s) Vt, s^O. 

In this way we may define, for any s ^ 0 , the measure Ug as the image of under the 
transformation x i—)• V(s)x 

i^s{dO = (V(s)#//t(s)) (dO V s ^ 0 , 
where ^ stands for the push-forward operation on measures. 


/ = [ 4>{V{s)x)iJ,t(s)idx) V(/> G Cb(M), s ^ 0. (3.20) 

</ M IR. 

Notice that, whenever nt is absolutely continuous with respect to the Lebesgue measure over 
M with utidx) = f{t,x)dx, the measure Ug is also absolutely continuous with respect to the 
Lebesgue measure over M with i/ 5 (d^) = 5 '(s,^)d.^ where 

9is,0 = '^’(■s)■^ f{t{s),V{s)~^() G M, s ^ 0, 

which is nothing but (1.6). Such a definition of Ug allows to define the semi-flow {J-g)g^o which 
given any initial datum /xq £ satisfying (1.14) associates 

= i^s = V{s)#St(^g){^o) Vs^O, 

where {St)t^o is the semi-flow associated to equation (1.1). The decay of the moments given by 
Theorem 3.5 readily translates into the following result. 

Theorem 3.7. Let 7 > 0 and let hq G 'P®(M) n 'Pexp, 7 (l^) be a given initial datum. Denote by 
iFg{fj,o) = Ug for any s ^ 0. Then, 


(1) When 7 G (0,1], there exists some universal eonstant JC^ > 0 (not depending on hq) sueh 
that for any p ^ 7 


min <j M..y (po)^ ; ( ^ ) }^Mp{ug) Vs^O, 


and for any p ^ 2, 

Mp {ug) if max <j Mp{iao ); 
(2) When 7 > 1, for any p > 1 


2p \ 

1-aP -bP 


Vs ^ 0. 


min <{ Ml {fiof ; ( - ^ j | ^ Mp (ug) V s ^ 0 , 

C-y,iA^ V ^ 


(3.21a) 


(3.21b) 


(3.22a) 



20 


R. ALONSO, V. RAGLAND, Y. CHENG & B. LODS 


and for any p ^ 2 

p 

Mp{us) Vs^O, (3.22b) 

where and Ai are defined in the above Proposition 3.4- 
Proof. The proof follows simply from the fact that 

Mp{us) = V{srMp{iit(^,)) ysfiO, 

where {pt)t is the weak measure solution to (1.1) associated to po. Then, according to Theorem 
3.5 and using that exp(ps) = (l + 'jt{s)Y^'^, we see that for 7 G (0,1] it holds 


f M^(^o)(l + 7 t(s)) \ 7 

Vl + IC.y M.y{po)t{s) J 


if Mp (lys) ^ Mp{po) 


1 + 7t(s) 

1 + ^{l-aP -bP) Mp {po)i t{s) 


P 

7 


where the lower bound is valid for any p 'J while the upper bound is valid for any pfi 2. Since 
min(l, ^) ^ ^ max(l, ^) for any A,B,t ^ 0 we get the conclusion. We proceed in the 

same way for 7 > 1 . □ 


An important consequence of the above decay is the following proposition. 
Proposition 3.8. Let po G ® given initial eondition. Assume that 

Mp{po) ^ Mp ypfi max( 7 , 2 ), 

where 

P 

2p 


(3.23) 


M„ , 

^ ' 1 - aP - 6P 


y p max( 7 , 2 ). 


Then, pQ G 'Pexp, 7 (l^) o-fid there exists an explieit a > 0 and C = C{a) > 0 sueh that 

sup [ exp(a|^P)i/^(dO ^ C, 

s70 JR 

where Ug = Fs{po) for any s 0. 


Proof Let us first prove that pQ G Vexp^'yi^)- Notice that for any 2 : > 0 

r p 

/ exp(2;|^|^)^o(dC) = • 

JR p^g P- 

Let us denote by po the integer such that 7 Po ^ max( 7 , 2) and 7(^0 ~ 1) < max( 7 , 2). Using 
Stirling formula together with the fact that limp_^oo(l ~ — b'^P) = 1, one can check that there 

exists some explicit a > 0 such that the series 

p=po 

converges for any 0 ^ 2 ^ a which gives the result. Now, we may define Ug = J-'s{po) for any 
s ^ 0. As previously, for any 2 > 0 



DISSIPATIVE BOLTZMANN EQUATION IN 1-D 


21 


and we deduce from Theorem 3.7 that 

Mp{us) ^ Mp Vp ^ max( 7 , 2), s ^ 0 . 

The conclusion follows. □ 


Remark 3.9. ITe do not need to derive the equation satisfied by Vg since we are interested only 
in the fixed point of the semiflow. However, using the fact that actually provides a strong 

solution to (1.1), using the chain rule it follows that 

[ [ 0(6/^o(dO+ / dr / + [ (Q (i^r , ; </>)dr Vs^O, 

Jm. m Jo Jm. Jo 

(3.24) 

for any (j) G (IR) and where <p' stands for the derivative of fi. 

The link between solution to (1.13) and the semiflow (^" 5)570 is established by the following 
lemma. 


Lemma 3.10. Any fixed point pL G 'P^(M) n 'Pexp 7 (l^) of the semi-flow {J-s)s^o is a solution to 
(1.13). 


Proof. Let /r be a fixed point of the semi-flow (Ts)s, that is, Hs{pi) = pL for any s ^ 0. Then, 
according to (3.20), for any G Cfe(M) 


[ <('(0/^(d0 = / fi{V{s)x)pt(s){dx) 

Jm, Jm. 

where = 5t(s)(/r). In particular, choosing s = s{t) 

[ = [ (j){V{t)x)ixt{dx) 

J M. J'M. 

Applying the above to 4>(V{t)~^ ■) instead of 4>, one obtains 

/ 4>{Vit)~^^)tJ-idf,) = [ 4>{x)ixt{dx) 
J'M. Jm. 


y s 0, 


ytfio. 


ytfio. 


Computing the derivative with respect to t and assuming (j) G C^(M), we get 

;^(ni)“') [ ^<P'{y{t)-^^)KdO = {Qil^t,Pt); fi) ytfiO. 
at 

Using the definition of V(t) it finally follows that 

-{l + 'yt)~^~^ [ {Vityas,) tx{dC) = {Q{pt,Pt)] (p) ytfiO. 

Jr 


Taking in particular t = 0 it follows that pi satisfies (1.13). 


□ 


4. Existence of a steady solution to the rescaled problem 

4.1. Steady measure solutions are steady state. In this section we prove Theorem 1.5, 
that is, we prove that steady measure solutions to (1.9) are in fact functions provided no mass 
concentration happens at the origin. The argument is based on the next two propositions. 

Proposition 4.1. Let pi G ^max(7 2)(^) “ steady solution to (1.9). Then, there exists H G 

L^(M) such that 


Cpi{df) = H{OdC. 
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Proof. Introduce the distribution (3{^) := ^/x(d^) which, of course, is defined by the identity 

/ /3(0V'(0dC = / for any V’ G • 

J M. Jm. 

One sees from (1.13) that /3 satisfies 

= Qifi, fj.) (4.1) 

in the sense of distributions. Since G ’P(^ax (7 2 )(^)’ follows that /i) belongs to 7W'’'(M). 

Therefore, as a solution to (4.1), the measure /3 is a distribution whose derivative belongs to 
A4(M). It follows from [19, Theorem 6.77] that (3 G BViocfM.) where BV(S.) denotes the space 
of functions with bounded variations. This implies that the measure (3 is absolutely continuous. 
In particular, there exists H G L^(M) such that /3(d^) = This proves the result. □ 

Proposition 4.2. Let /r G ^max (7 2 )(^) ® steady solution to (1.9). Then, there exist kq ^ 0 

and G G L^^_^ 2 ) (1^) nonnegative sueh that 

= G{f)df + Ko (Jo(dO 

where 5o is the Dirae mass in 0. 


Proof. Let us denote by B{M.) the set of Borel subsets of M. According to Lebesgue decomposition 
Theorem [31, Theorem 8.1.3] there exists G G L^(M) nonnegative and a measure fig such that 

fi{do = G{m+f^sm 

where the measure fig is singular to the Lebesgue measure over M. More specifically, there is 
T G B{M.) with zero Lebesgue measure such that fig(S.\T) = 0. The proof of the Lemma consists 
then in proving that fig is supported in {0}, i.e. T = {0}. This comes directly from Proposition 
4.1. Indeed, by uniqueness of the Lebesgue decomposition, one has 


ep(de) = eG(c)d^ + ep,(do = H{f)df, 

with PI G L^(M), so that 

C/^s(dO = 0. 

This implies that /r^(M\(—5; 5)) = 0 for any (5 > 0 and therefore that fig is supported in {0}. 
Notice that, since fl G ^max( 7 , 2 )(^) ^ ^ -^max( 7 , 2 )(^)- 

Proof of Theorem 1.5. With the notations of Proposition 4.2, our aim is to show that kq = 0. 
Plugging the decomposition obtained in Proposition 4.2 in the weak formulation (1.13), we get 

- /e0'(OG(Odc= [ Q(G,G)(c)</>(Ode + No / icr(</>(aO + </>(&e)-</>(e)-</>(o))G(e)de, 

Jm Jm Jm ^ ' 

where we used that Q((Ioj<^o) = 0 = / ^i;^^(0'^o(d0- Recall the hypothesis 

Jr 


^7= / ierG(e)de>o, 
Jr 


from which the above can be reformulated as 

Koin^m= [ QiG,G){mm 

Jr 

+ / ^./-'(^^(OdC + Avo / |^r(<(>(aO + </>(H)-</>(6)G'(6de 
Jr Jr ^ ^ 


(4.2) 
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for any (j) G Notice that the above identity can be rewritten as 

Jm 

for some -functions 


^ G 


Nn 

+ -r 


a 

\aj 

b 

b \ 


and -B(^) = Let (phe a smooth function with support in (—1,1) and satisfying 0(0) = 1. 

For any e > 0, 0 (j) belongs to one can apply (4.3) to get 

No = / ryl(^)0 G(C)-0'd^. 


Hence, 

0 < No ^ ||0||z,oo /" |yl(O|d^-hsup|C0'(OI / G(0d^- 

J—£ —e 

Letting e —)■ 0, one obtains no = 0, thus, using hypothesis (4.2) we must have no = 0. □ 

4.2. Proof of Theorem 1.6. We have all the previous machinery at hand to prove the existence 
of “physical” solutions to (1.9) in the sense of Definition 1.3 employing the dynamic fixed point 
Theorem 1.1. Let us distinguish here two cases: 

1) First, assume that 7 G (0,1]. Setting y to be the space A4(M) endowed with the weak-* 
topology, we introduce the nonempty closed convex set 

Z:=||UGP°(M) such that / |?r/u(dO ^ , and f ^ Mp Vp ^ 2| 

I Jr ^7 Jr j 

where Mp was defined in Proposition 3.8 and JC^ > 0 is the positive constant given in Theorem 
3.7. This set is a compact subset of y thanks to the uniform moment estimates (recall that y is 
endowed with the weak-* topology). Moreover, according to Proposition 3.8, there exists a > 0 
and C{a) > 0 such that, for any £ Z, 

[ exp(Q;|,f0) /u(d^) ^ C{a) < 00 . 

Jr 

Thus, fi G Vexp,'r{^)- Therefore, using Theorem 1.4, {St{fi))t^o and, consequently, (J^s(/i))s^o 
are well-defined. Setting Ug = BgilJ'), it follows from Theorem 3.7 that 


[ |C|Vs(dC)^Mp Vp^2, s^O, 

Jr 


Using the lower bound in (3.18a), we deduce that 

[ |Cr*^s(dC) ^ ^ Vs ^ 0. 

Jr ^7 

This shows that Ug £ Z for any s ^ 0, i.e. J'siZ) C Z for all s ^ 0. Moreover, one deduces 
directly from Proposition 2.9 that (V^s)s is continuous over Z. As a consequence, it is possible to 
apply Theorem 1.1 to deduce the existence of a measure fx £ Z such that a steady 

measure solution to (1.9) in the sense of Definition 1.3. Finally, since fx £ Z, its moment of order 
7 is bounded away from zero and by Theorem 1.5, // is absolutely continuous with respect to 
the Lebesgue measure. This proves the result in the case 7 G (0,1]. 
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2) Assume now 7 > 1 and let T = max( 7 , 2). Then, consider y to be the space A4(M) endowed 
with the weak-* topology and we introduce the nonempty closed convex set 




such that f |,f|^ |u(d^) ^ Mp Vp ^ T, and f |C|/u(d^) ^ £ 
^ M. M 


for some positive constant i to be determined. In fact, we prove that there exists i = £( 7 ) 
sufficiently small such that Z is invariant under the semi-flow {J-s)s^o- Indeed, according to 
(3.22a), for any £ > 0 if po is such that Mi(/io) ^ then 


Mi(JT(po)) ^ mini i; 


^7,1 A " 


Vs > 0. 


where > 0 is some positive universal constant. And, according to (3.9), Ai is any real 
number larger than max / C-f^i ; where is another universal positive constant. 


number larger than max < 1 ; -;—r-r;— >, w 

1 Mi(po)i+^i 

In particular, choosing i small enough such that 


it is possible to pick Ai := , where we recall that, since hq ^ Z and 1 -|- 7 ^ max{ 7 , 2}, 

one has Mi+^(po) ^ In such a case, one gets 




= min < £ ; 


We set I in order to get 


min < i ; 




I J 

and Mi[Ts{lJ'o)) ^ ^ for any s ^ 0. Arguing as in the case 7 G (0,1], this shows that Ts{Z') C Z 
for any s ^ 0 and, there exists a steady measure /i which is absolutely continuous with respect 
to the Lebesgue measure. 

5. Numerical simulations 

This section contains numerical simulations for the rescaled equation 

dsg{.s,i)-\{d^ s^O, (5.1) 

where Q{g,g) has been previously defined in ( 1 . 2 ) and (1.3). We recall that such a model has 
been studied in [17, 29] in the case of 7 = 0 and it admits a unique steady state 


such that 




Mi({) e d{= 0 
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We will use the numerical solutions of (5.1) to verify the properties of our models for general 
values of 7 . We shall consider here initial datum go{0 — which shares the same first 

moments of Mi i.e. 





1 

0 

1 


(5.3) 


The coefficient c = — ^ (a^+ 6 ^ —l) = aft > 0 in equation (5.1) is the only one that gives stationary 
self-similar profile with finite energy in the case 7 = 0, see [17]. In contrast, as already noticed, 
the case 7 > 0 accepts any arbitrary positive coefficient in the equation, thus, we will perform all 
numerical simulations with such coefficient for comparison purposes. Recall that in our previous 
analysis we choose this coefficient to be 1 . 


5.1. Numerical scheme. To compute the solution, we have to make a technical assumption 
which is the truncation of M into a finite domain = [—L, L], When L is chosen large enough 
so that g is machine zero at ^ = ±L, this will not affect the quality of the solutions^. Then, we 
use a discrete mesh consisting of N cells as follows: 

-L = ^1 < ^3 < ... < = L . 

We denote cell Ij = (Cj_i jCj+i)) with cell center and length A^j = 

^ ., 1 — The scheme we use is the discontinuous Galerkin (DG) method [18], which has 

•^2 ^2 

excellent conservation properties. The DG schemes employ the approximation space defined by 

= {Vh : Vh\i, G PHIj), 1 ^ J ^ N}, 

where P^{Ij) denotes all polynomials of degree at most k on Ij, and look for the numerical 
solution gii G , such that 


[ dsgh{s,C)vh{0<iC + h{a^+ b‘^^9h{s, () d^VhiO 

+ \{a^ + h^- 1 ) ({ighvl)j_i - {ighv^)-^i^ 

= [ Qigh,gh){s,()vh{0<i^, 

JI 3 

holds true for any G V^. In (5.4), ^g^ is the upwind numerical flux 


j = l,...,N (5.4) 


/ ?%(s,0 if? 

I ?5'ft(s,?) if ?<0 


where g^,g^ denote the left and right limits of g^ at the cell interface. Equation (5.4) is in fact 
an ordinary differential equation for the coefficients of gh{s, ?). The system can then be solved by 
a standard ODE integrator, and in this paper we use the third order TVD-Runge-Kutta methods 
[30] to evolve this method-of-lines ODE. Notice that the implementation of the collision term in 
(5.4) is done by recalling (1.3), and we only need to calculate it for all the basis functions in V^. 
This is done before the time evolution starts to save computational cost. 

The DG method described above when k ^ 1 (i.e. we use a scheme with at least piecewise 
linear polynomial space) will preserve mass and momentum up to discretization error from the 


^Although we have no rigorous proof that the tails of g will decay, at least, as fast as Ae , this is expected. 
For dimensions d > 1, this is a well known fact for the elastic and inelastic Boltzmann equations. 
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boundary and numerical quadratures. This can be easily verified by using appropriate test 
functions Vh in (5.4). For example, if we take u/j = 1 for any j, and sum up on j, we obtain 

f dsOh = T - l) {9h {i = L)+ gl = -L)) . 

J 

If L is taken large enough so that gh achieves machine zero at ±L, this implies mass conservation. 
Similarly, we can prove 

[ dsgh^dC =-^{a^+ b‘^-1) [ s-h+ ^(«^ + &^ - 1) (s-h (? =-^^) - =-i)) • 

Jn Jn ^ 

Again, when L is large enough and the initial momentum is zero, this shows conservation of 
momentum for the numerical solution. 

5.2. Discussion of numerical results. We use as initial state the discontinuous initial prohle 


5(0,6 


if|?l<v^ 

0 otherwise. 


This prohle clearly satishes the moment conditions (5.3). We take the domain to be [—40,40] 
and use piecewise quadratic polynomials on a uniform mesh of size 2000. Four sets of numerical 
results have been computed, corresponding to ( 7 , 0 ) = (1,0.1), (1,0.3), (1,0.5), (2,0.5), and 
(3,0.5) respectively. The computation is stopped when the residual 


/ f f 9h{s'^+\0-gh{s^,0 

Vn[ As 

reduces to a threshold below 10 “^ indicating convergence to a steady state. 

In Figure 1 we plot the objects of study in this document, that is, the equilibrium solutions 
for different values of 7 . In this plot, the amplitude of the solutions has been normalized to one 
at the origin for comparison purposes. The numerical solutions are used for the cases 7 = 1,2,3, 
while for 7 = 0, we use the theoretical equilibrium Mi as dehned in (5.2). In general terms, these 
smooth patterns are expected with exponential tails happening for any 7 > 0. The behavior of 
the prohles at the origin is quite subtle and will depend non linearly on the potential, for instance, 
the case 7 = 2 renders a wider prohle relative to 7 = 0 in contrast to 7 = 1 or 7 = 3. This is 
not to say that such behavior is discontinuous with respect to 7 , it is simply the net result of 
the contributions of short and long range interactions of the particles in equilibrium. 

In Figure 2, we hx 7 = 1 and compare the stationary solution for different values of a. 
Recall that the parameter a measures the “inelasticity” degree of the system with a = 0 being 
elastic particles and with a = 0.5 being sticky particles. As expected, smaller values of a will 
render a wider distribution prohle at the origin keeping the tails unchanged. Near the origin, 
the distribution of particles for less inelastic systems will be underpopulated relative to more 
inelastic systems which force particles to a more concentrated state. Tails, however, are more 
dependent to the growth of the potential and should remain relatively unchanged despite changes 
in inelasticity. Interestingly, the numerical simulation shows a unexpected effect: the maximum 
density of particles is not necessarily located at the origin. 

In Figure 3, we plot the evolution of energy as a function of time in a system of sticky 
particles a = 0.5 using different values of 7 . Changes in the relaxation times are expected since 
the potential growth 7 impacts directly on the spectral gap of the linearized interaction operator. 
This numerical result seems to conhrm, in one dimension, the natural idea that higher 7 implies 
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Figure 1. Rescaled equilibrium solutions for different values of 7 with a = 0.5 
(sticky particles). Curves corresponding to 7 = 1,2,3 are computed numerically, 
while the curve for 7 = 0 is obtained from the known steady state Mi in (5.2). 





Figure 2. Equilibrium solutions for different values of inelasticity a when 7 = 1. 

higher spectral gap, hence, faster relaxation to equilibrium. Refer to [28] for ample discussion in 
higher dimensions for the so called quasi-elastic regime. Additionally, the results of Figure 3 are 
numerical confirmation of the optimal cooling rate given in our Theorems 3.5 and 3.7. 

In Figure 4, we investigate the evolution of the distribution function and its discontinuities for 
the case 7 = 1 and a = 0.5. The simulation shows, in our one dimension setting, a well stablished 
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Figure 3. Evolution of energy for 7 = 1 , 2 , 3, when a = 0.5. 

phenomena happening in elastic and quasi-elastic Boltzmann equation in higher dimensions: 
discontinuities are damped at exponential rate [28]. As a consequence, points of low regularity 
which are contributed by Q^(g, g) due to such discontinuities will be smoothed out exponentially 
fast as well. This is the case for the point ^ = 0 in this particular simulation. A numerical 
simulation was also performed using an initial Gaussian profile (not included). Both numerical 
simulations showed an evolution towards the same equilibrium profile which confort us in the 
belief that the constructed self-similar profile is unique. 

Finally, we verify the performance of our scheme by plotting the distribution’s mass, momen¬ 
tum and energy. Only the plot for 7 = l,a = 0.5 is shown since the other cases display similar 
accuracy. In Figure 5 we plot the evolution of mass, momentum, energy and residual (in the log 
scale). The decay of residual shows convergence to steady state, while mass and momentum are 
preserved up to 10 digits of accuracy verifying the performance of the DG method. 

6. Conclusion and perspectives 

In the present paper we studied the large time behavior of the solution to the dissipative 
Boltzmann equation in one dimension. The main achievement of the document is threefold: (1) 
give a proof for the well-posedness of such problem in the measure setting, ( 2 ) provide a careful 
study of the moments, including the optimal rate of convergence of solutions towards the Dirac 
mass at 0 in Wasserstein metric, and (3) prove the existence of “physical” steady solutions in 
the self-similar variables, that is, steady measure solutions that are in fact absolutely continuous 
with respect to the Lebesgue measure. Let us make a few comments about the perspectives and 
related open problems. 

6.1. Regularity propagation for inelastic Boltzmann in 1-D. The numerical simulations 
performed in Section 5 seem to confirm that many of the known results given for inelastic Boltz¬ 
mann in higher dimensions should extend to inelastic Boltzmann in 1-D, at least, under suitable 
conditions. More specifically, rigorous results about propagation of Lebesgue and Sobolev norms, 









ene 
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Figure 4. Evolution of g{s, when 7 = 1,0 = 0.5. 




Evolution of energy Evolution of residual 




Figure 5. 7 = 1, o = 0.5, evolution of mass, momentum, energy and residual. 
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and exponential attenuation of discontinuities for the time evolution problem should hold. Sim¬ 
ilarly, the study of optimal regularity for the stationary problem is an interesting aspect of the 
equation which is unknown. 

6.2. Alternative approach d la Fournier-Laurengot. Exploiting the analogy between (1.8) 
and the self-similar Smoluchowski’s equation, one may wonder if the approach performed by N. 
Fournier & Ph. Laurengot in [21] can be adapted to (1.9). We recall that the approach in [21] 
consists in finding a suitable discrete approximation of the steady problem for which a discrete 
steady solution can be constructed. If such discrete solution exhibits all the desired properties 
(positivity, uniform upper bounds and suitable lower bounds) uniformly with respect to the 
discretization parameter, then, one can pass to the limit to obtain the desired steady solution to 
(1.9). Such approach fully exploits the 1-D feature of the problem. Besides, it does not resort 
to the evolution equation (1.8), fact that makes it very elegant. The main contrast with respect 
to [21] lies in the fact that no estimates for moments of negative order seem available for our 
problem. Moreover, Smoluchowski’s equation is such that the collision-like operator sends mass 
to infinity while the drift term brings it back to zero. The model (1.9) has the opposite behavior: 
the collision tends to concentrate mass in zero while the drift term sends it to infinity. 


6.3. Uniqueness and stability of the self-similar profile. Now that the existence of a steady 
solution to (1.9) has been settled, the next challenge is to prove that such self-similar profile is 
unique and that it attracts solutions to (1.8) as s —)• oo or, at least, to find conditions for this to 
hold. This is certainly the case in the simulations performed in Section 5 which show, in addition, 
exponential rate of attraction. For the 3-D inelastic Boltzmann equation, such a result has been 
proven in [28] in the so-called weakly inelastic regime (a perturbation of the elastic problem). 
Since the 1-D Boltzmann equation is meaningless for elastic interactions a perturbative approach 
seems inadequate. Once a uniqueness theory is at hand, it would be desirable to obtain rate of 
convergence, see for instance [3]. This would render a more complete picture of the large time 
behavior of the dissipative Boltzmann equation on the line. 


6.4. The rod alignment problem by Aranson and Tsimring. Aranson and Tsimring in [4] 
have introduced the following model for rod alignment (the rods have distinguishable beginning 
and end) 


dtP{t,e) 



+ |0*|^d0*-p(t,0) j" p{t,9 +e,)\e,yde, ( 6 . 1 ) 


having initial condition P{0) = Pq, angle 6 £ (—vr, tt). The authors introduced the model for 
Maxwellian interactions 7 = 0 , yet, the model is sound for any 7^0. We refer to [8, 16] for 
other variations of such model. Here P{t,9) is the time-distribution of rods having orientation 
9 £ [—7r,7r). Equation (6.1) models a system of many discrete rods aligning by the pairwise 
irreversible law 

[9 — ^,9 + ^ (^ 9 , 9 ). (6.2) 

Let us explain the interaction law (6.2). We start by fixing a horizontal frame and picking two 
interacting rods with orientation 0i,02 £ [~7r, tt). Define 0* G [— tt, vr) as the angle between the 
ends of the rods. Bisect the rods and define 9 £ [—vr, vr) as the angle between the horizontal 
frame and the bisecting line. Thus, we can express the rods orientation, up to modulo 27r, by 
the relation 9i = 9 — ^ and 92 = 9 + ^ with respect to the horizontal frame. After interaction, 
both rods will align with the bisection angle 9. This law produces the alignment of rods, we refer 
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to [4, 8 ] for an interesting discussion and simulations. The law (6.2) can be written in terms of 
the rod orientations 9i, 62 ^ [—vr, vr) as 


{01,02) 


( ( 61+62 6i+02'J 

\9i- 92 

[ (^+7r’^+7r) 

\9i-02 


^ vr, 
> vr. 


(6.3) 


Note that in the case |0i — ^ 2 ! >71 the addition of tt is needed since we chose the alignment 
to occur in the direction of the bisecting angle associated to the ends of the rods (as opposed 
to the beginnings of the rods). The interaction law (6.3) is discontinuous, thus, intuitively we 
understand that model ( 6 . 1 ) will not have conservation of momentum because there is a choice 
of alignment direction. Let us fix this by considering an initial datum Pq with compact support 
in (—7r/2,7r/2) 

SuppPo C (- 7 r/ 2 , 7 r/ 2 ). 

Such a property is conserved by the dynamic of (6.1) and it corresponds to a system of rods where 
rod’s beginning and end are indistinguishable, thus, we can always assign an angle 9 £ (— 7 r/ 2 , 7 r/ 2 ) 
to each rod. For such a model, the weak-formulation is very similar to that of (1.1) except for 
the fact that all integrals are considered now over the finite interval (— 7 r/ 2 , 7 r/ 2 ). For this reason, 
the decay of the moments of the solution P{t) to (6.1) is identical to that of (1.1). Consequently, 
this translates into the convergence of P{t) towards a Dirac mass centered at 0 as t —)• 00 in the 
Wasserstein metric. The question is to understand the model after self-similar rescaling where 
the support of solutions is no longer fixed and given by ( — V{t)TT/2,V{t)7r/2) —)• (— 00 , 00 ) as 
t — 7 - 00 . Thus, it is natural to expect that the self-similar solution to (6.1) will converge towards 
the steady solution to (1.9). 


6.5. Extension to other collision-like problems. It seems that the present approach is ro¬ 
bust enough to be applied to various contexts. In particular, the argument may be helpful to 
tackle notoriously difficult questions, such as, the existence of a stationary self-similar solution 
to the Smoluchowski equation with ballistie kernel interaetions. It may be possible, also, to give 
a more natural treatment of the stationary inelastic Boltzmann equation in the framework of 
probability measures. The difficulty will be to find a dynamical stable set in order to apply the 
dynamical fixed point and a suitable regularization theory for the stationary equation of the 
particular problem. 


Appendix A. Cauchy theory in both the L^-context and the measure setting 


In this Appendix, we give a detailed proof of the existence and stability estimates of Section 
2.2 yielding to Theorem 1.4. We fix here 7 > 0 and set T = max( 7 , 2 ). We begin with an 
existence and uniqueness result for the Cauchy problem (1.1) in the special case in which the 
initial datum is absolutely continuous with respect to the Lebesgue measure, i.e. 

;Uo(dx) = /o(x)dx. 

Theorem A.l. Fix (5 > 0. Let a nonnegative fo £ (M) be given with /o ^ 0. Setting 

/io(dx) = /o(x)dx, there exists a unique family {ft)t^o C L 2 _,_.^^ 5 (M) sueh that /if(dx) = ft{x)dx 
is a weak measure solution to (1.1) assoeiated to hq. Moreover, 


/t(x)dx = / /o(x)dx, / xft{x)dx = / x/o(x)dx, 

^ Jm. Jm. 

[ |xp/t(x)dx ^ [ |xp/o(x)dx Vt^O. 


and 


(A.l) 
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In addition to this, if one assumes that 


foe[]Ll(R), 

then {ft)t^o C 

Proof. We follow here the approach of some unpublished notes of Bressan [15]. For K > 0 and 
5 > 0, we introduce 

= jo ^ / G f f{x)dx= I fo{x)dx, I /(x)|xp+'^+^dx ^ kI . 

For / G ^k,5 ^ changing variables in the collision operator leads to 

IIQ(/>/)IIli ^ [ [ fix)f{y)\x-y\'^ {l + {bx + ayf)dxdy 

JR JR 

+ [ [ f{x)f{y)\x-y\'^ {l + x‘^)dxdy. 

</ M. M. 

Now, for any {x,y) G 

|x - yP ^ C^{\xP + lyD 

with C.y = max{l; Therefore, one observes that 

II Qif, f)hi ^ M^(/)Mo(/) + 2C^ M 2 (/) M^if) + 2C^ M 2 +^(/)Mo(/), 

where M^if) = f^f{x)\x\^dx for any k ^ 0. Notice that Mo{f) = ||/o|| li. Using Young’s 
inequality, for any 0</c<2 + 7 + (5 one has Mi ^{ f ) ^ K + ||/o||li from which we deduce that 
there exists C > 0 such that 

\\QifJ)\\Ll^C{\\fo\\lr+K^) . 

Consequently, Q(/,/) G L 2 (M). Let us now prove that the restriction to of the mapping 
/ I—S' Qif, f) is Holder continuous. For f,g^ klK,s, 

l|Q(/,/) - Q(ff,5)llLi ^ / / \if - g)ix)\fiy)\x-y\'^il + ibx + ayf)dxdy 

JR JR 

+ [ [ 9ix)\if-g)iy)\\x-y\'^il + {bx + ayf)dxdy 


+ [ [ \if - 9 )ix)\fiy)\x-y\'* {l + x‘^)dxdy 

</ M. */ IR. 

+ [ [ 9ix)\if - g)iy)\\x-y\^ {l + x^) dxdy . 

M. */ IR. 

Proceeding as previously, one notices that there exists C > 0 such that 

l|Q(/, /) - Qi9,9)\\Ll ^ CiWfoiy + K)\\f - ghr + 2C'^||/ohi / |(/ - g)ix)\\xf+^dx . 

J M 

Thanks to the Holder inequality, we have 

[ \if - 9)ix)\\x\‘^^'^dx if ( f |(/-y)(x)||xp+^+'^dxy (f |(/- y)(x)||xpdxy 

»/]R \t/IR / \»/]R / 

^(2iL)W^||/-y||2f . 

^2 
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Combining the previous two inequalities, we deduce that the mapping / i—)■ Q{f, f) is uniformly 
Holder continuous on when restricted to Let us look for a one-sided Lipschitz 

condition. For f,g ^ we introduce 


[/, g_ _ 


lim 

5 ^ 0 ” 


11/ + 'Sfl'ILi “ II/IIli 

s 


The dominated convergence theorem implies that 

[/,5]_ ^ [ sign{f{x))g{x){l + x‘^)dx. 

Jr 

Our aim is to show that there exists a constant L > 0 such that for any f,g£ ^k,Si 


[f - g.QifJ) - Q{9,9)]_ ^ L\\f - 9\\li- 


But, 


/ sign((/- 5 r)(x))(Q(/,/) - Q{g,g)){x){l + x‘^)dx 
Jr 

\{f- 9 )ix-ay)\{f + g){x + by)\y\^{l + x‘^)dxdy 

^ Jr Jr 

+ A / [ {f + g)ix-ay)\{f-g){x + by)\\y\'^{l + x‘^)dxdy 

^ Jr Jr 

+ a// \{f- 9 ){x + y)\{f + g)ix)\y\'^{l + x‘^)dxdy 

^ Jr Jr 

-III \if - 9 ){x)\{f+ g){x + y)\y\'^{l + x^)dxdy. 

^ Jr Jr 

Thus, changing variables leads to 
[/ - 5, Q(/,/) - Q{9,9)]_ 

^11 I \if - 9 )i^)\if + 9 ){y)\x - y\'^ {3 +{ax + byf + {bx +ayf+ y^) dxdy 

^ Jr Jr 

-A / / \{f- 9 ){x)\{f + g){y)\x-yp{l + x^)dxdy. 

^ Jr Jr 

Finally, we obtain 

[/-5, Q(/,/) - Q(ff,5)]- ^ I I \if-9){x)\{f + 9){y)\x-y\^{l + y'^)dxdy 

Jm. Jm. 

^ mf -9 \\li > 

with L = 2C'..),(4||/o||/,i -|- 3K). Next, let us look for a sub-tangent condition. Given / G ^k,s 
and h ^ 0, one notices that 


f{x) + hQ{f,f){x) = h [ f{x - ay)f{x+ by)\y\'^dy + f{x) (l-h[ f {x + y)\y\'^ dy 

In particular, what prevents / -|- hQ{f, f) to be a.e. nonnegative is the influence of large x in the 
last convolution integral. To overcome this difficulty, for any R > 0, we introduce the truncation 
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fnix) = f{x)x{\x\<R}- Then, since f fn one deduces from the above identity that 
f{x) + hQifn, fR){x) ^ h [ fnix - ay)fR{x + by)\y\'^dy 


+ fR{x) 1 - /i / /r(x + y)\y\^dy a.e. x € 


(A.2) 


Now, 


fR{x + y)\y\'^dy = [ //j(y)|x - ypdy < max( 2 '^ \l) [ f{y){\x\'^ + \y\'^)dy 
J'M. Jm. 


and using Young’s inequality, one sees that there exists some positive constant Cq > 0 depending 
only on K,5,x and ||/o||li but not on R such that 


/ fR{x + y)\ypdy ^ Co{l + |xp) for any x G M, i? > 0 and / G ^k^s- 
Jr 

Therefore, recalling that /r is supported on {|x| < R}, one deduces from (A.2) that 

fix) + hQ{fR,fR){x) ^ 0 for a.e. X G M, \f0 < h < Hr := ^ . 

Moreover, since Q preserves the mass, 

/ {fix) + hQifR,fR)ix))dx= f /(x)dx= f /o(x)dx. 

Jm. Jm. Jm. 

Finally, using (2.7) it follows that, for any i? > 0 

[ QifR,fR)ix)\x\'^"^^^^dx 
Jr 

^(l - 0^+^+'^ - 6'^+^+^) f [ /R(x)/ij(y)|x-yp^+^+'^dxdy ^ 0. 

^ ^ Jr Jr 

Consequently, 

/ (/(x) +/iQ(/ij,/ij)(x))|xp+^+‘^dx ^ [ /(x)|xp+^+^dx ^ a: VA > 0 . 

</ M. </ M 

We have thus shown that, for any R > 0 and any 0 < h < Hr, one has / + hQ{fR, /r) G ^k,s- 
In particular, for any R > 0 and any 0 < h < Hr one has 

dist(/ + hQ{f,f),nK,5) 

^ 11/ + hQif, /)-(/ + hQifR, fR)hi = h II Qif, f) - QUr, fR)hi . 


Now, for / G ^k,5i one can make ||Q(/,/) — Q(/r,/r)||li arbitrarily small provided i? > 0 is 
large enough so that the sub-tangent condition 

liminf h~^dist (/ -h hQ(/, /), Q,k,s) = 0 
h—>- 0 + 

holds true. We may now apply [26, Theorem VI.4.3] and deduce the existence and the uniqueness 
of a global solution / to (1.1) such that /(t) G i^K,5 for every t ^ 0. Moreover, (A.l) holds and, 
if /o G Tfc(K), it follows from (2.7) that for every k ^ T and every t ^ 0, 
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which implies (together with the conservation of the mass) that f{t) G for any t ^ 0 and 

any /c G M. Finally, it is easily checked that the family dehned by nt{dx) = f{t,x)dx for 

any t ^ 0 is a weak measure solution to ( 1 . 1 ). □ 


Proof of Proposition 2.7. The proof follows the approach of [25, Section 4] and we only sketch 
the main steps of the proof. First, since /Uq ^ Dirac mass centered at 0, the 

temperature Tq := JjgX^;Uo(dx) is positive and one can dehne a sequence (Fq^)„) such that 


lim [ (p{x)FQ{x)dx = [ (p{x)fiQ{dx) V y? G L“x(^) D C(M). (A.3) 

Jr Jr 


with Fq G PI^^qLJ(M) for any n ^ 1 (notice that, as in [25], Fq is some slight modihcation of 
the Mehler transform of //q)- Then, according to Theorem A.l, for any n ^ 0, there exists a 
family (T”)f^o C Ly(M) such that (^f)t^o is a weak measure solution to (1.1) associated to ^q, 
where 


Then, noticing that 
one easily checks that 


p.f{dx) = Fj^{x)dx Vt^O. 

^ iii^nix ^ iimoIIt vn^i, 


; p )\ ^ 2||v9||oo||iioIIt Vt ^ 0, n ^ 1. 

from which one deduces as in [25] that there exists C = C{p.Q) > 0 (depending only on H/toIIt) 
such that, for any t 2 ^ ^ 0 , 


sup 

n^l 



p{x)pJl^{dx) 


p{x)PtMx 


^ C'(iUo) ||<^||oo 1^2 - tl| Vv 2 GCb(M). 


(A.4) 


Moreover, on the basis of the a priori estimates (2.10) (see also Remark 2.5), 

MkiFn^CkijAmonpUrt-"^ yk>2. 

Since lim„_^oo ||/^ollo = 1 and lim„_^oo lll^ollx = ||/^o||x according to (A.3), we deduce that, for 
any A: > 2, there exists some positive constant Ck depending only on k, 7 , and Mx(/io) such 
that 

sup Mk{Ff') ^ Ckt~^ V/c>2. 

n^l 

From this, we conclude, as in [25] that there exists a subsequence (still denoted by) and 

a family (^j )t70 c Afx(i^) such that 


lim 

n^oo 


p{x)Pt(.<Jx) 


(p{x)iJ.t{dx) V(/2GCc(M), t^O, 


(A.5) 


ll/Uillx ^ IImoIIx , Mfc(//t) ^ Cfct r yt> 0 , k> 2 , 

and (A.4) still holds for the limit fj,t (which implies that, for any ip G Cfe(M), the mapping 
t G [ 0 ,00) I—)• Jjg (/2(x)/ri(dx) is continuous). To prove that {iit)tpo is a measure weak solution 
to (1.1) associated to pQ in the sense of Dehnition 1.2, one argues exactly as in [25, Section 
4]. Finally, the fact that hq G nfc 70 ^fc(^) implies that {iit)tpo C nfe 70 ^fc(^) proven as in 
Theorem A.l. □ 


Remark A.2. Arguing as in [25], it is not difficult to prove that any weak measure solution to 
(1.1) is in faet a strong solution in the sense of [25]. 
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Proof of Proposition 2.8. Let T > 0 be fixed. For any (j) G Lip]^(M) and t G [0, T] define 

:= f (j){x)iJ.tidx) - f (j){x)ut{dx). 

IR. M 

We will also use in the proof the notation 


W(t) = dKRipt, vt) = sup W^{t ), 

0eLipj(R) 


and recall that 


W{t)= \x - y\TTt{dx,dy) for some m G ■ 

Jr'2 

Let now f G Lip^(M) be fixed. One has 
d 1 /* 

= o / I® - vVMix, y)yt{dx)yt{dy) 
dt ^ 2 J^2 

— - / \v — w\^/S.(l>{v,w)vt{dv)vt{dw) 

2 yR2 

= [ 7rt(dx,dr;) [ \\x - y\'^Ao(f{x,y) - \v - w\'^Ao(f{v,w)\TTt{dy,dw) 


(A. 6 ) 


(A.7) 


(A. 8 ) 


where Ao(/>(x,y) = (j){ax + by) — (j){x) for any {x,y) G M^. Now, 

\x — yPAo(?i>(x, y) — \v — w\'^Ao(j){v, w) = 

- yP - lu - w\'^'^Ao(j){x,y) + |u - w\'^ (^Ao(p{x,y) - Ao(f{v,w)^ , 

and, recalling that the Lipschitz constant of cj) is at most one, the second term readily yields 

AQ(f{x,y) — Ao(j){v,w) = (j){ax + by) — (j){av + bw) — (p{x) + 4 i{v 

^ (1 + a)|x — ul + b\y — w\ . 

For the first term in (A. 7 ),we use the identity A = min{A, B} + (A — B)^, valid for any A,B ^ 0 , 
to obtain the estimate 

\A"< - B'^\A =\A^ -B"i\[ min{A, B] + {A- B)+) 

^ 1^7 _ B"i\ min{A, B] + |A^ - B'^\ |A - B\ 

^ (1 + 7) max{A, B}'^\A — B\. 

The last inequality follows noticing that \A"< — B'^\ min{A, i?} ^ 7 max{A, A — B\ . Since 
|Ao(/)(x, y)\ ^ b\x — y\ we can choose A = \x — y\ and B = \v — w\ to conclude that 


(^\x — y\'^ — \ v — Aq4>{x, y) ^ 6(1 + 7 ) max{|x — y\,\v — w\}'^\\x — y\ — \ v — w\ 

^ 6(1 + 7 ) max{|x — y\,\v — rc|}"’'(|x — v\ + \y — rc|) . 


(A.9) 


Gathering the estimates (A.7),(A. 8 ) and (A.9) in (A. 6 ) and using symmetry of the expression, 
it follows that ^ + 7)77(t) where we introduced 

H{t) := [ 7 rt(dx,du) [ (Ix^ + lyP + |x - u| 7 rt(dy, drc). 

Jr‘2 Jr2 


(A.IO) 
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Expand H(t) = Hi{t) + H 2 {t) where 


Notice that 


Hi{t) : = [ TTt{dx,dv) f [\y\'^ + \w\'^)\x - v\TTt{dy,dw) 

Jr2 Jm.2 

H2{t) : = [ (|xP + |i;P )\x - p| 7 rt(dx, dv). 

Jr‘2 

Hi{t) = J^[\y\'^ + \w\'^) 7 rt{dy,dw) y Jx - p| 7 rt(dx,dp) 

= W{t) [ \x\'^{fit + i^t){dx) ^CW{t). 

Jr 


(A.ll) 


The last inequality follows because the weak measure solutions yt and ut have the 7 -moment 
uniformly bounded in t G [0)2^]) and additionally, VTt G achieves the Kantorovich- 

Rubinstein distance. We estimate now H 2 {t) as in [22, Corollary 2.3]. Namely, for any t G [0,T] 
and any r > 0 , one has 


/ ( |xp-|-[pp ) |x — p| 7rt(dx, dp) ^ 2r'^ / |x — p| 7rt(dx, dp) 

JR^ \ ) Jr^ 


+ / ( |xP-|-| pP j |x — p| 7 ri(dx, dp) 

= 2 r'''lT(t)-|- f Mx|'>'-|-[pp^x — p| 7 ri(dx, dp), 

since irt G ut) achieves the Kantorovich-Rubinstein distance. Setting now such that 

(|xP-h|pp) |x-p| +expr^i|L^^ 


^ Re (exp(e|xP) -|- exp(e|p|'>')) V (x, p) G 


it follows that 


Choosing 


we obtain 


H 2 it) ^ 2r'^W{t) + Re Ct{£) exp ( 


r'*' = | 21 oglT(t)/e| 


H 2 {t) ^ -W{t)\logW{t)\+ReCTi£)W{t) 


(A.12) 


Estimates (A.ll) and (A.12) imply that 


-W^t) < KeCrie) W{t){l + [log W(t)|) , 


(A.13) 


with a constant Ke > 0 depending only on 7 and e > 0. Integrating (A.13) and taking the 
supremum over cj) G Lip^(M) we get the conclusion. □ 
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